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We study some general aspects of the formalism of potential NRQCD (pNRQCD), an effective field theory that 
deals with ultrasoft degrees of freedom in Heavy Quarkonium systems. Specific attention is paid to its effective 
Lagrangian that it is displayed at the present level of accuracy. 
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1. Introduction 

The last years have witnessed important 
progress in the theoretical understanding of 
heavy-quark-antiquark systems near threshold 
through the use of effective field theories [y-Q 
(see also §-|l§ for related work in QED). The 
key point relies on the fact that, since the quark 
velocity t; is a small quantity, w ^ 1, a hierarchy 
of widely separated scales: m ^ mv S> mv^ ... is 
produced in these systems, where m is the heavy 
quark mass (hard scale), mv the soft scale and 
mv'^ the ultrasoft scale. One can take advan- 
tage of this hierarchy by integrating out the scales 
above the energies we want our theory to describe, 
i.e. above the ultrasoft ones. 

After integrating out the hard scale Non- 
relativistic QCD (NRQCD) is obtained 0. The 
Lagrangian of NRQCD can be organised in pow- 
ers of I/to. After integrating out the soft scale 
in NRQCD, potential NRQCD (pNRQCD) is ob- 
tained P]. The Lagrangian of pNRQCD is or- 
ganised in powers of 1/m and the relative coor- 
dinate r (multipole expansion). The matching 
coefficients in this latter case are non-analytical 
functions of r. 

The integration of degrees of freedom is done 
in practice through a matching procedure (see 
|[0,|,01 for details). The matching from QCD to 
NRQCD can always be done perturbatively since, 
by definition of heavy quark, m 3> Aqcd @i|lli ■ 
The matching from NRQCD to pNRQCD can 
only be carried out perturbatively when mv ^ 
^QCD- We will assume this to be so throughout 
this work. Therefore, the matching coefficients in 
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both NRQCD and pNRQCD can be computed or- 
der by order in as . The non-analytical behaviour 
in 1/to appears through logs in the matching co- 
efficients of the NRQCD Lagrangian: 

771 

c~ Ausiln h-B), 

where /Jh denotes the matching scale between 
QCD and NRQCD. A typical matching coefficient 
for pNRQCD could have the following structure: 

D - V{r, p, si, S2) {A' \nmr + B'lnfir + C), (I) 

where /j, denotes the matching scale between 
NRQCD and pNRQCD and V denotes a non- 
analytic function of r. 

2. Theoretical Framework 

2.1. pNRQCD: the degrees of freedom 

Integrating out the soft scale, mv, in NRQCD 
produces pNRQCD [Q. The relevant degrees of 
freedom of pNRQCD will depend in general on 
the nonperturbative features of NRQCD. In this 
work we will assume that there exists a matching 
scale /x such that m,v ^ /z ^ m,v^, Aqcd, where 
a perturbative picture still holds. 

Strictly speaking pNRQCD has two ultravio- 
let (UV) cut-offs Ai and A2. The former fulfils 
the relation mv^, Aqcd ^ Ai <^ mv and is the 
cut-off of the energy of the quarks and of the en- 
ergy and the momentum of the gluons (it corre- 
sponds to the fj, above), whereas the latter fulffis 
m,v ^ A2 ^ m and is the cut-off of the rela- 
tive momentum of the quark-antiquark system, 
p. In principle, we have some freedom in choos- 
ing the relative importance between Ai and A2. 



Our choice is K\/m <Si Ai, which guarantees that 
the UV behaviour of the quark propagator in pN- 
RQCD is that of the static one. 

If we denote any scale below Ai, i.e. 
ttiv'^jAqcd, •■■, with Amp, we are in a posi- 
tion to enumerate the effective degrees of free- 
dom of pNRQCD. These are: Q-Q states with 
energy of 0{Amp) and relative momentum not 
larger than the soft scale, and gluons with en- 
ergy and momentum of 0{Amp)- Let us define 
the centre-of-mass coordinate of the Q-Q system 
R = (xi + X2)/2 and the relative coordinate 
r = xi — X2. A Q-Q state can be decomposed 
into a singlet state S'(R, r, t) and an octet state 
0(R, r,t), in relation to colour gauge transfor- 
mation with respect to the centre-of-mass coordi- 
nate. We notice that in QED only the analogous 
to the singlet appears. The gauge fields are eval- 
uated in R and t, i.e. A^^ = ^/^(R, t): they do not 
depend on r. This is due to the fact that, since 
the typical size r is the inverse of the soft scale, 
gluon fields are multipole expanded with respect 
to this variable. 

2.2. Povirer counting 

In order to discuss the general structure of the 
pNRQCD Lagrangian, let us consider in more de- 
tail the different scales involved in the problem 
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The variables r and m will explicitly appear in 
the pNRQCD Lagrangian, since they correspond 
to scales that have been integrated out. Although 
p and 1/r are of the same size in the physical 
system, we will keep them as independent. This 
will facilitate the counting rules used to build the 
most general Lagrangian. 

With the above objects, the following small di- 
mensionless quantities will appear: 
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(2) 



m rm 

Note that l/pr or Amp/p are not allowed since 
p has to appear in an analytic way in the La- 
grangian (this scale has not been integrated out). 
The last inequality in Eq. (||) tells us that r can 
be considered to be small with respect to the re- 
maining dynamical lengths in the system. As a 



consequence the gluon fields can be systemati- 
cally expanded in r (multipole expansion). 

Therefore, the pNRQCD Lagrangian can be 
written as an expansion in 1/m (from the first two 
inequalities of Eq. (^), and as an expansion in r 
(the so-called multipole expansion, corresponding 
to the third inequality of Eq. (||)). 

2.3. pNRQCD: the Lagrangian 

The pNRQCD Lagrangian would look as fol- 
lows: 

-CpNRQCD = Tr -( St ( 00 

+Ot(^zi?o-^-V;(°)(r) + ...jO| 

-hgVA(r)Tr{0V-ES + SV-E0} 
+g^^^^TT {O^r • E O + OtQr • E} 
_^pa^p,.a^ (3) 

where the dots indicate higher-order potentials in 
the 1/m expansion, we have neglected center-of- 
mass variables and we have just kept 0(r) terms 
in the multipole expansion. Let us now display 
the structure of the potentials up to 0(l/m^). 
Order l/?7i°. By dimensions V} (r) can only 
have the following structure: 
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Order 1/m. By dimensions plus time reversal 
Vs (r) can only have the following structure: 

l^._££i£^, c. = «. (5) 

m Zmr^ 

Order l/ni^. The following matching potentials 
appear to the accuracy to which the matching has 
been performed at present 
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Note that p appears analytically in the poten- 
tials. The power on n (p") to which p appears in 
the potential is constrained by the power in \/va. 
Some of the operators above become too singular 
once the dependence on In r is taken into account 
in the matching coefficients. For instance, one 
could have 5{y) In mr which is ill defined. In fact, 
this operator has to be understood as the sum of 
two operators, well defined each: 5{y) In ^^ and 
reg\ + In/i/j. The latter is basically the Fourier 
transform of In ^ (see |^ and references therein). 
The matching potentials for the octet sector 
could be defined in a similar way by changing the 
overall factor C/ by C/ — Ca/2 and the subscript 
s by o in each of the matching coefficients (ay^ , 

The above representation of the potential can 
be related with others found in the literature [|2| 
by the use of the equations of motion. 

In order to have the proper free-field normal- 
isation in the colour space we define S = -j^S 

and O = ^0°, where Tp = 1/2, and S and O 
are 1/2 x 1/2 tensors in spin space. 

In principle, terms with higher time derivatives 
acting on the fields 5' or O could be considered 
in the pNRQCD Lagrangian. These terms are re- 
dundant in the sense that one could make them 
disappear from the Lagrangian and still correctly 
predict all physical observables by re-shuffiing the 
values of the matching coefficients. We have 
chosen the minimal form of the pNRQCD La- 
grangian, where higher time derivatives are ab- 
sent. In fact, one can always get rid of these 
terms by systematically using local field redefini- 
tions without changing the physical observables 
(spectrum) of the theory. 

2.4. pNRQCD matching coefficients 

In order to obtain the different matching coeffi- 
cients in pNRQCD: ay,, D'^'^\ D^^) ... one has to 
perform the matching between NRQCD and pN- 
RQCD. A detailed description of the procedure 



can be found in P,|9|,|7|,[l3| . Here we just state the 
main results and briefly discuss how they are ob- 
tained. The matching coefficients read as follows 
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{r)cl ^ as{r) (l - Ca^ In mr) .(14) 

In order to obtain ay^ with the accuracy above, 
it is necessary to perform the matching between 
NRQCD and pNRQCD (at 0(l/m°)) exactly at 



the two-loop level and with the leading-log accu- 
racy at the three-loop level, i.e. to compute the 
static potential to this order. The one-loop re- 
sult was obtained in Ref. Q, the two-loop one 
in Ref. Q and the three- loop leading-log in ||]. 
The Pn are the coefficients of the beta function 
and the values of ai and 02 can be found in Ref. 
]15| (see also [|l5[ for notation). 

For ni ' the first contribution appears at one 
loop and, to our knowledge, was first calculated 
in pffl . The two- loop leading-log term is obtained 

in |l3| (see also |17 ). 

(2) 
The tree-level result for the different Ds can 

be obtained from the results of positronium (see 
H). The one- loop leading- log dependence on 
In fir is obtained in [Of (see also |17|). By taking 
into account the NRQCD matching coefficients 
(c£), cp ...; see [p|,pl| for its values at one loop) 
in the vertices when matching to pNRQCD, one 
obtains the leading- log dependence on In mr plus 
also some finite pieces. In any case, since the full 
0{al) correction to the different D's is not yet 
available, the results above should only be trusted 
up to finite pieces at one loop. It is also worth 
noting that the spin-dependent matching poten- 
tials do not suffer from ultrasoft divergences at 
one- loop. 

In order to obtain the leading-log dependence 
on n of the matching coefficients above it is 
enough to compute the UV divergences of pN- 
RQCD at next-to-leading order in the multipole 
expansion and up to 0(l/m^) Q (see also [p7[). 

With the matching coefficients above, the 
leading-log, 0{mal), corrections to the heavy 



quarkonium spectrum can be computed |13 17[ |. 

With little effort a few things can be said with 
respect to the octet matching potentials, ay^ co- 
incides with av^ at one loop (although not nec- 
essarily beyond^) and also the leading-log (three- 
loop) dependence is known M . We are not aware 
of any calculation of Do (of which the first non- 
zero contribution also appears at one- loop). The 

(2) 
tree-level contribution to the different Do can 

also be read for the positronium calculation, keep- 
ing in mind that, unlike in the singlet case, in 
this case the annihilation diagram gives a nonzero 



contribution (as in QED). 

The leading-log dependence of Va and Vb on /x 
is also known [0. It reads 

8 a 
VA{r,fi} = VB{r,fi) = 1 + -CA—lnrii. (15) 
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